Three-Band Model for Quantum Hall and Spin Hall Effects 
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Topological properties of a certain class of spinless three-band Hamiltonians are shown to be 
summed up by the Skyrmion number in momentum space, analogous to the case of two-band 
Hamiltonian. Topological tight-binding Hamiltonian on a Kagome lattice is analyzed with this 
view. When such a Hamiltonian is "folded", the two bands with opposite Chern numbers merge 
into a degenerate band exhibiting non-Abelian gauge connection. Conserved pseudo-spin current 
operator can be constructed in this case and used to compute the pseudo-spin Hall conductance. Our 
model Hamiltonians belong to the symmetry class D and AI according to the ten-fold classification 
scheme. 

PACS numbers: 



I. INTRODUCTION 

Topological properties embedded within a band struc- 
ture have emerged as one of the central themes of con- 
densed matter physics nowadays. The topological nature 
is often expressed as the gauge field derived from the un- 
derlying wave functions and manifests itself in such ob- 
servable phenomena as anomalous and spin Hall effects 
in metals 1 ' 2 . Quantized versions of the phenomena ex- 
ist in insulators, when a certain topological number can 
be associated with the completely filled bands. Classifi- 
cation of the permissible topological numbers in a given 
band Hamiltonian was recently carried out 3 . Numerous 
efforts are being made at the moment to suggest specific 
models belonging to a particular entry in the classifica- 
tion table with relevance to realizable experimental sys- 
tems. A most striking instance of this kind of effort is 
Haldane's proposal of microscopic two-band Hamiltonian 
with spontaneous quantized Hall conductance 4 , its gen- 
eralization to time-reversal-invariant quantum spin Hall 
phase 5,6 , and its subsequent realization in HgTc/CdTc 
heterostructure 6,7 . 

Basically, spontaneous quantum Hall models and a 
subset of quantum spin Hall models defined as two 
time-reversal copies of the quantum Hall model can be 
phrased as a property of the generic two-band Hamilto- 
nian H = V'k^kV'k, where ipk is a two-component 
spinor and the matrix Ji^ is Jik = e^l + dk • er (<x=Pauli 
matrix). As long as the dk-vector remains nonzero 
over the entire Brillouin zone (BZ) such that a unique 
unit-vector c?k = dk/|dk| exists, the two eigenstatcs 
of 'Hk will be classified according to their spin helicity 
dk ■ O'lipu) = ^iV-'k )■ From the Abclian gauge connection 

ak = -£(V>k l^klV'k) = +*(^k l^klV^k ) tne associated 
flux density can be written in two equivalent forms, 



Pk 



da y 
dk T 



daa_ 

dk,, 



1 t / <9d k dd k \ 



Quantized number for the band, obtained as the integral 
(l/27r) J dk x dk y of Eq. (f), can therefore be interpreted 
as either the first Chern number or the Skyrmion number 
of the given band 8-10 depending on the use of the second 
or the third term of Eq. (f ) for the integrand. 

There are quite a few other models embodying non- 
trivial topological numbers that involve the use of 4x4 
r-matrices, many of which are summarized in Ref. 3. Al- 
though the general classification scheme makes no ref- 
erence to the dimensionality of the matrix itself, in 
practice almost all explicit examples of band Hamilto- 
nians with nontrivial topology take the form of an even- 
dimensional matrix in momentum space. In sharp con- 
trast, general discussion of the topological character for 
odd-dimensional 3x3 matrices, or three-band Hamiltoni- 
ans, appear to be lacking. The work of Ohgushi, Mu- 
rakami, and Nagaosa (OMN) 11 provided, to the authors' 
knowledge, the first microscopic example of a three-band 
model Hamiltonian with non-trivial topological Chern 
number. The three-sublattice structure of the Kagome 
lattice is a natural platform for the three-dimensional 
Hamiltonian matrix to arise. 

In this paper, we divide the general three band Hamil- 
tonian into the spin-1 part and nematic part then try 
to solve each part separately. The spin- 1 model and ne- 
matic model involve quantum Hall effect and spin Hall 
effect respectively. In contrast to two band model, for 
the spin-f model, the first Chern number is one half of 
skyrmion number. We realize that the existence of half- 
skyrmion (mcron) configuration for a unit Chern number 
band. When we consider the Kagome lattice models for 
the nematic model, which the spin matrices represent the 
pseudo-spin from the three-sublattice structure, our spin 
Hall effect describes the pseudo-spin Hall effect. 

The remaining part of this paper is organized as fol- 
lows. In Sec. II, in terms of spin-I matrices, we consider a 
particular three band model involving the quantum Hall 
effect. All of the other three band model is topologically 
equivalent to spin-1 model unless the band gap is closed. 
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As a realistic example, we introduce the OMN model on 
Kagome lattice which has unit Chern number including 
the meron structure in momentum space. In Sec. Ill, by 
introducing the nematic operator, we construct the dif- 
ferent type of three band model with time-reversal sym- 
metry which involves the (pseudo-)spin Hall effect. For 
an specific model based on Kagome lattice, we perform 
the linear response calculation of the (pseudo-)spin Hall 
conductivity. Then we give an physical interpretation 
for our (pseudo-)spin Hall effect. In Sec. IV we conclude 
with brief summary and discussions. 

II. SPIN-1 MODEL 

In general, arbitrary three-band Hamiltonian can be 
written as a linear combination of eight Gell-Mann ma- 
trices A a (a = l,-- ,8), H k = £kI + ELi d k' Aa > charac- 
terized by eight-component real field a\. Due to its com- 
plicate structure, the exact analytical solution of general 
three band model is difficult to obtain. By the way, as is 
well-known, Gell-Mann matrices are the generators of the 
SU(3) Lie algebra which possess as subgroups several sets 
of three matrices forming SU(2), or SO (3) sub-algebra. 
In particular the subset that generates the SO (3) rotation 
are the matrices of spin-1 operators. We will show that 
the subset of 3x3 Hamiltonians spanned by the SO(3) 
generators can be readily analyzed in analogous fashion 
as their 2x2 counterparts. We specialize to three-band 
Hamiltonians of the type 

H k = d k • S, (2) 

where one possible choice of the spin-1 matrix S would 
be {S a ) M = -ie af3y , or S = (A 7 , —A 5 , A 2 ): 

/o o o\ /o o A 

S x = A 7 = -i , S y = -A 5 = 0, 
\0 i J \-i 0/ 

/0 -i 0\ 

S z = A 2 = \i . (3) 
\0 0/ 

All other choices of SO (3) subgroup elements ought to 
be related to this by a suitable unitary transformation. 
As in the two-band Hamiltonian, we assume |dk| remains 
nonzero throughout the BZ. The eigenstates of energies 
dk and can be worked out readily in terms of the unit 
vector dk = (d x , d y , d z ) 12 ' 13 , 

^ / d x d z + id y \ / d x \ 

h£)= , = d y d z -id x , K)= d v -( 4 ) 

V2[i-(4) 2 ] V dl-l, ) \dj 

The state of energy — dk is the complex conjugate of | ip£ ) : 

h£> = (W£>)*- 

With the aid of explicit wave functions one can evalu- 
ate the gauge flux associated with each band. The cen- 
tral band, having real- valued wave functions, has the zero 
flux, while the upper band has the flux density 12 



(da y da x \ ~ ( dd k <9d k \ 

p * = {w x -wJ = dk - [dk- x x w y )- (5) 

The lower band has the opposite sign of the flux density. 
Comparing Eq. (5) to Eq. (1) one notes a difference of 
factor 2, originating from spin-1 being twice the size of 
spin-1/2 12 . A full Skyrmion for dk in the BZ will thus 
imply an even Chern number for the spin-1 bands. 

We ask now if certain topological three-band models 
considered in earlier literature can be framed in the form 
of the "parent Hamiltonian" dk • S plus some perturba- 
tion. As long as the parent model is connected smoothly 
to the full one without the gap closing the topological 
property will be completely captured by the relation, Eq. 
(5). The Kagome lattice model of OMN 11 is 

4 OMN) =d k -S + d^.S', (6) 

where S = (A 7 , -A 5 , A 2 ), S' = (A 6 , A 4 , A 1 ), while d k = 
2sin[0/3](cos(k • a 2 ),cos(k • a 3 ),cos(k • ai)) and d^ = 
2cos[(/>/3](cos(k-a2),cos(k-a 3 ),cos(k-ai)) for the three 
unit vectors a, (i = 1, 2, 3) of the Kagome lattice shown in 
Fig. 1(d). The flux <f> penetrating the triangle unit of the 
Kagome lattice is compensated for by —20 flux through 
the hexagon in the OMN model. We have verified that 
d^ • S' does not close the energy gap of the parent Hamil- 
tonian dk • S, provided the flux (f> satisfies the condition 
sgn (sin[(/>/3]) = — sgn(sin(/)). The Skyrmion number as- 
sociated with the dk-vector of the parent OMN Hamilto- 
nian is readily obtained, N s = — (l/2)sgn (sin[</>/3]), giv- 
ing rise to the Chern number for the upper band C + = 
—sgn (sin[0/3]) according to Eq. (5). The Chern number 
obtained for the same band in Rcf. 11 is C + = sgn(sin (/)), 
equal to our result provided sgn (sin[0/3]) = — sgn(sin</>). 
This is precisely the same condition required for the topo- 
logical equivalence of the original OMN model to its par- 
ent form. Written in real space lattice the parent Hamil- 
tonian gives the hopping amplitude proportional to +i 
for every nearest neighbor bond traversed in the counter- 
clockwise direction. 

An astute reader may wonder how come the Skyrmion 
number in the parent OMN Hamiltonian is only half of 
an integer. In fact any Hamiltonian matrix dk • S that 
shares the periodicity of the first BZ, i.e. dk+G = dk 
for reciprocal lattice vector G, will only allow the inte- 
ger Skyrmion numbers due to the homotopy ^(S* 2 ) = Z 
and, by virtue of Eq. (5), only even Chern integers. On 
the other hand, explicit calculation of the Berry phase 
flux yields the Chern numbers of ±l n . So why the ap- 
parent contradiction? A closer examination shows that 
dk for the OMN model is not periodic under k — > k + G 
but rather under k — > k + 2G as shown graphically in 
Fig. 1(a). The same relation between Chern number 
and winding number for three orbital model with the 
spin-1 representation is studied in Rcf. 12. As a result, 
dk+G • S is equivalent to dk ■ S only up to some constant, 
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FIG. 1: (color online) (a) Skyrmion density (background 
color) for the parent Hamiltonian of the OMN model. Color 
bar represents the Skyrmion density. The dk-vector (arrows) 
is periodic over four BZs as indicated by the black hexagon. 
White hexagon is the BZ. (b) Energy dispersion for open 
boundary condition of the parent Hamiltonian dk ■ S. The 
one-dimensional BZ is indicated by two dashed vertical bars 
at k x = ±7r/2. Edge modes are shown in red, bulk modes in 
black, (c) Energy dispersion of the folded OMN Hamiltonian 
— (dk ■ S) 2 under the open boundary condition, (d) Open ge- 
ometry used in the calculation of (b) and (c). Hopping occurs 
along the solid bonds only. Three unit vectors ai,a2,a3 are 
shown. Three sublattices sites are indicated as A, B, C. (e) 
Spin Hall conductance (in units of e/2n) for the folded OMN 
model as a function of the Fermi energy. 



k-independent unitary rotation: [d k +G'S] = {/' [dk-S]Z7. 
By evading the periodic condition the Hamiltonian be- 
comes exempt from the usual homotopy consideration as 
well, making the half Skyrmion configuration possible. In 
fact with a different choice of the basis the OMN matrix 
can be written in such a way that "Hk+G = "Hk holds. In 
that case, however, the model can no longer be reduced to 
the form dk • S and the Skyrmion number interpretation 
of the Chern number fails to apply 14 . 

Topological nature of the spin-1 Hamiltonian %k = 
dk • S can be phrased in the general classification scheme 
of Ref. 3. Under complex conjugation one has S* = — S 



time- reversal symmetry (TRS) and/or particle-hole sym- 
metry (PHS) if some unitary transformation could re- 
late U%. to +H-k (TRS) or -"H-k (PHS), respectively. 
When the dk-vector is even, dk = d_k, we further have 
"H-k = "Hk, hence the TRS/PHS amounts to the ex- 
istence of unitary transformation U^1~L\JJ = T~Hk, re- 
spectively It follows that SLS is trivially satisfied with 
U = I, while TRS cannot be achieved with any U. This 
places our Hamiltonian (2) in the D class 3 , where the al- 
lowed topological numbers are the integers in two spatial 
dimensions, equal to the Chern numbers we just calcu- 
lated. If instead we had d k = dk the Hamiltonian 
would have TRS but not PHS, placing it in class AI, with- 
out any topological numbers. Indeed one can easily show 
that for the dk-vector of odd symmetry the Skyrmion 
number vanishes identically. 



III. NEMATIC MODEL 

So far the discussion of the topological character of 
three-band Hamiltonian is restricted to the "spin" type 
given by Eq. (2). The nomenclature is obviously derived 
from S being a representation of S — 1 spin. On the other 
hand, the following observation prompts us to study an- 
other class of three-band Hamiltonians that we denote 
the "nematic" type. Note that the anti-commutators of 
the three spin operators (S x , S v , S z ) generate the follow- 
ing five: 



S x S y + S V S X = A 1 



S V S Z + S z S y = A e 



S z S x + S x S z 



[s x ] 2 - [s y ] 



y}2 



v/3 



(2[S Z } 2 - [S x ] 2 - [S v ] 2 ) = A 




and %l = (d k ■ S)* 



-Hk- The Hamiltonian possesses 



These are precisely the remaining five of the Gcll-Mann 
matrices, or in the language of spin liquids, the ne- 
matic operators. An arbitrary 3-band Hamiltonian is 
therefore a sum of the spin part and the nematic part, 
%k = dk-S + Dk-N, where the five components of N refer 
to the above five operators and Dk is a five-component 
function of k. The Hamiltonian consisting solely of the 
nematic operators Dk-N are real, [Dk-N]* = Dk-N, and 
TRS/PHS conditions become C/t[D k • N]U = ±D_ k ■ N. 
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In this case even (odd) Dk = +(— )D_k generates the AI 
(D) class Hamiltonian 3 . 

A particularly simple kind of nematic Hamiltonian 
arises by folding the previous spin Hamiltonian 



•H k = d k -S^-H 



(/). 



[d k -S] 5 



(8) 



The folding ensures the degeneracy of Kramers' pairs 
l^k) w ith the identical energy — (dk) 2 wherein non- 
Abelian gauge connection arises 15-18 . The same wave 
functions worked out before the folding remain as eigen- 
states of the new Hamiltonian. A similar three-band 
Hamiltonian, with d k oc k, was proposed for the p-orbital 
bands of Si 18 . An example of the folded four-band model 
employing the 5 = 3/2 spin operator S is the Luttinger 
Hamiltonian, extensively studied by Murakami, Nagaosa, 
and Zhang as a model for dissipationless spin Hall cur- 
rent in GaAs 16,17 . Even for spinless models such as ours, 
an analogue of spin Hall current can be defined and its 
response function computed. The special case of dk oc k 
was analyzed in Rcf. 18 . Here we maintain the framework 
as general as possible by keeping dk an arbitrary non- 
zero vector over the BZ. 

Each band has the associated helicity number, dk • S = 
h (h = +1,0, —1) and the band-dependent flux density, 



F: 



h di. 



cWk 9c?k 



(9) 



The sum of Chern numbers for the doubly degenerate 
band is obviously zero. Regarding S as the pseudo-spin 
matrix, the projected pseudo-spin density may be defined 



as 



16-16 



~ a 



p 

,(0) 



<5 P 



J2 Pi a) s a p( a \ (10) 



with P p uj = 1 - (dp ■ S) 2 and P P L > = (d p ■ S) 2 denot- 
ing the projection onto the eigenstates with Eq = 
and Ei = — (d p ) 2 , respectively, and ip p is the three- 
component spinor consisting of A, B, C sublattice site 
operators. Together with the projected spin current op- 
erator 



s a , 



dn P f) 

dp 



(11) 



they obey the continuity equation = iq ■ Jq + 
0(q 2 ). The d.c. pseudo-spin Hall conductivity follows 
as (V^volumc) 17,18 



AO) 



(12) 



where is the Fermi function of each band a = 0,1. 
Here j stands for the direction of the applied electric 
field, i the spatial direction of spin current, and a is the 
spin orientation. This formula is general and applicable 
to any three-band Hamiltonians H£ of the folded form, 
Eq. (8). Numerical evaluation of the pseudo-spin Hall 
conductance a% y for the parent OMN model is shown in 
Fig. 1(e). Due to symmetry other spin orientations give 
the same Hall conductance: a xy = a xy — u xy 

Intuitive understanding of the pseudo-spin Hall con- 
ductance follows from the fact that spin operator we use 
also serves as the sublattice current operator. For in- 
stance, S x expressed in the sublattice basis |A), |B), |C) 
becomes S x = — i|B)(C| + i|C)(B|, equal to orbital cur- 
rent operator Jg C for that bond. As a whole we may 
identify S = (•/bc^cA'^ab) m obvious notation. The 
eigenstate |^) has the average ("0k l^lV'k) = ^k, and 
with our new interpretation it implies that the current 
loop around a ABC triangle shown in Fig. 2(a) is 

(^Wab + Jbc + JcaM£) = d x + d y + d z and mi- 
nus this value for \tp k ). These opposite current loops are 
the analogues of opposite spin orientations in the genuine 
spin Hall effect. Now these current loops of opposite signs 
move in the opposite directions due to the opposite signs 
of the gauge flux they experience, conserving the time- 
reversal invariance, and result in the pseudo-spin Hall 
phenomena (Fig. 2(b)). 




FIG. 2: (color online) (a) Two opposite current loops around 
a ABC triangle, (b) Pseudo-spin Hall phenomena. Two op- 
posite loops move in the opposite directions. The electric field 
is applied along vertical axis. 

Imposing boundaries in the y-direction at y = and 
at y = L y as shown in Fig. 1(d) introduces some 
edge modes. The topologically protected edge modes 19 
of the parent OMN model connecting the topological 
bands (top and bottom) with the non-topological band 
(center) as shown in Fig. 1(c) obey the dispersion 20 
Ek x = ±2 cos k x when the hopping magnitude is chosen 
to one. The folded Hamiltonian — (dk ■ S) 2 also support 
edge modes, albeit unprotected in the topological sense, 
with the dispersion E^ x = — (3 + cos2fc a; ) as shown in 
Fig. 1(d). The explicit derivation of the edge dispersion 
is presented in the Appendix. For the OMN model the 
sign of the velocity dE^ / dk x determines the spatial lo- 
cation of the edge as cither y ~ or y ~ L y . In the 
folded case both edge modes obey the same energy dis- 
persion. The absence of protected edge mode connecting 
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the degenerate bands to the flat band is due to the net 
topological number for the former bands being zero 19 ' 21 . 

IV. CONCLUSION AND DISCUSSION 

In this paper we have taken up the study of the topo- 
logical aspects of three-band Hamiltonian, applicable to 
spinless fcrmions and bosons with three-fold orbital or 
sublatticc degree of freedom. As with the topological 
two-band models, the Hall conductances of the three- 
band Hamiltonian is shown to be governed by the same 
(Chern number) <-> (Skyrmion number) correspondence 
with a critical difference of factor two allowing for the ex- 
istence of meron structure in momentum space for three- 
band Hamiltonians. For the OMN model we succeeded in 
obtaining the topological number of the model by using 
the parent model of Eq. (2). However, we didn't mention 
about how to obtain the parent model from the general 
three band model. In this point, we need further studies. 
A second class of three-band Hamiltonians obtained by 
the folding procedure is studied. The double degeneracy 
ensured by folding naturally leads to non-Abelian gauge 
structure and spin Hall phenomena. Both of three-band 
models we study, before and after the folding, are subject 
to the symmetry classification scheme of Ref. 3. 
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Appendix: Edge modes in the folded model 

In the OMN model, the chiral edge state is obtained 
in Ref. 20. Here we calculate the edge mode solutions of 



the folded OMN model, 

^ /} =-[d k -S] 2 , (A.l) 

where 

dk = 2 sin[0/3](cos(k • a2), cos(k ■ a^), cos(k • ai)). 

(A.2) 

The real space expression of the Hamiltonian is obtained 
by writing down its Fourier transformation, 

H {S) = - ^ ( <4[a r+ 2a 3 +Q r+ 2ai] 
r 

+ b\ [br+23.! + kr+2a 2 ] + 4 [ c r+2a 2 + C r +2a 3 ] 
Ct r [& r _|_ a3 _ a2 -f- &r+ai H~ Cr+a3 ^r+ai — a.2 ] 
L?r [c r _|_ a2 Cr+ag— ai H~ ^r+as— a 2 ^r+ai ] 

- 4 [a r +ai -a 2 + a r +a 3 + ^r+a 2 + 6 r +a 3 -ai ] J + h.C. 

(A.3) 

In order to obtain the edge states, we should choose a 
particular boundary condition. The boundary condition 
is depicted in Fig. 1 (d). Since the system is no longer 
periodic along the y direction, k y is not good quantum 
number. Let us consider momentum representation in x 
direction 



a(r) - a ni ,„ 2 = -= £ e- ik ^ +n ^a n2 (k). (A.4) 



Here, we used r = (2ni + ri2, v3ri2)- Inserting the one 
particle state 



i*(*)> = E (^w4w i w +^ c -|Wi|W) io> (a.5) 

j=l 

into the Schrodingcr equation H\H>) = E\^), we have the one-dimensional chain equations of ipj 19 

E^ = 2 cos k [i: a 3+ i + r 3 -i ~ e ik tf+L - e~*N#_$ - e ~^+i ~ e ^^-\) > ( A - 6 ) 
Etf.i =2cos(2fc)Vti +e lfe ^a + e^'VLi 

■J 2 •> 2 J 2 J ' 2 

-2e-* fe cosfcV?-i ~ 2e5 fc cosfcV? -2cosfcVLi -^ c -s., ( A -<0 

M_! =2cos(2fc)^_i +e- lk ^j c . 3 + e lk i' c i+1 

J 2 J 2 J 2 ■» 2 

- 2ei fc cosfcV?-i ~ 2e-^ k cosk^ - 2coskib b -x -ip h ,,i -ip h ,_ 3 , (A.8) 

J 2 3\ 2 J 2 
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with the boundary conditions 

2 2 2 

ib*+,=ib b ^=^b c ,i=0. (A.9) 

r max+l v max +i v max +^ v y 

For exponentially localized solutions on the boundaries 
we take the ansatz 22 

r 3 = rt l p, = »rV, i - »rV , 

(A.10) 

with the condition I?/) 6 ! = |?/> c |. Thus we can write as 

<f> b = l, ^ c = e* Vb,c = ve ieb >°, (A.ll) 
I 




where 77 is real function of momentum fc. Putting this all 
together we obtain the edge solutions as 

J J 2 ■* 2 

(A.12) 

with 

77= -(cosfc) ±1 , -Bfc = -(3 + cos2fc). (A.13) 

The degenerate edge spectrum are depicted in Fig. 1 (c) . 
For 7] 2 > l(rf < 1), the edge state is localized at up- 
per (lower) edge. 



FIG. 3: (color online) Plot of rf . For rf > l(r) 2 < 1), the 
edge state is localized at upper(lower) edge (blue(red) line). 
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